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ABSTRACT 

Aims. The circumstellar environment of L 2 Pup, an oxygen-rich semiregular variable, was observed to understand the evolution of 
mass loss and the shaping of ejecta in the late stages of stellar evolution. 

Methods. High-angular resolution observations from a single 8 m telescope were obtained using aperture masking in the near-infrared 
(1.64, 2.30 and 3.74 pm) on the NACO/VLT, both in imaging and polarimetric modes. 

Results. The aperture-masking images of Li Pup at 2.30 pm show a resolved structure that resembles a toroidal structure with a major 
axis of ~140 milliarcseconds (mas) and an east-west orientation. Two clumps can be seen on either side of the star, ~65 mas from the 
star, beyond the edge of the circumstellar envelope (estimated diameter is ~27 mas), while a faint, hook-like structure appear toward 
the northeast. The patterns are visible both in the imaging and polarimetric mode, although the latter was only used to measure the 
total intensity (Stokes I). The overall shape of the structure is similar at the 3.74 pm pseudo-continuum (dust emission), where the 
clumps appear to be embedded within a dark, dusty lane. The faint, hook-like patterns are also seen at this wavelength, extending 
northeast and southwest with the central, dark lane being an apparent axis of symmetry. We interpret the structure as a circumstellar 
torus with inner radius of 4.2 au. With a rotation velocity of 10 km s _1 as suggested by the SiO maser profile, we estimate a stellar 
mass of 0.7 M Q . 

Key words, stellar evolution - AGB stars - high-angular resolution techniques - Li Pup 


1. Introduction 


During the late stages of evolution of intermediate-mass stars 
(1-8 M G ) - that is, at the transition from the asymptotic gi¬ 
ant branch (AGB) to the post-AGB phase - the ejected mate¬ 
rial is shaped into intriguing forms (e.g., bipolar nebulae, torii, 
disks and/or spirals^ that depart from the initial spherical symme¬ 
try (Ivan Winckelll2003l and references therein). Although there 
are many studies in the literature of asymmetries found in post- 
AGB stars, detections of such asymmetries in AGB stars are rare. 
It is therefore uncertain at which point the asymmetries begin to 
develop. 

One of the brightest AGB stars in the near-infrared 
(iiik < -2.1 mag) near the solar neighborhood (64 ± 4 pc; 
Ivan Leeuwenl 2007) is L 2 Pup (HD56096, HIP34922, IRAS 


07120-4433). It is a semiregular variable (M5III; P = 140 + 1 
days) that experienced a dimming event starting around 1995 
(AV ~3 mag; Fig. lA.lk which was attributed to circumstellar 


dust dBedding et al.ll2002i ). The decline lasted for six years, and 
as of June 2014, the star has not yet begun to recover its former 
brightness. Dimmin g eve nts are present in a few other evolved 
stars (e.g., V Hya; iKnann et al.l 11999 ). which may repeat over 
the course of many years. L 2 Pup had a less severe dimming 
event around 1960; recovery took until 1980. Unlike V Hya, 
there is no clear periodicity to the dimming in the data, which 
go back to 1930. 

Only_a few imagesof L 2 Pup are available in the litera¬ 
ture. iBedding et al. ( 2002 ) stated that the IRAS images are non- 
descriptive. The WISE satellite observed L 2 Pup post-dimming 
in the near- and mid-infrared at low angular resolutions, although 
all images were saturated by the stellar emission. IJura et al.l 
(1 2002 ) studied L 2 Pup with the Keck 10 m telescope in the mid- 
infrared and inferred that two asymmetries are present at po¬ 
sition angles 135° and 225° east-of-north. Past studies on the 
variation of polarization of L 2 Pup are indicati ve of an asym- 
metricQ dust distribution (Magal haes et al.l fl986). In addition to 


* Based on observations from ESO programs 090.D-0541 & 090.D- 
0677 


1 In this paper, the term ‘asymmetric' will describe any aspherical 

structure. 
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the stellar component, llreland et al. ( 200413 ) inferred from the 
MAPPIT/AAT (3.9 m) observations in the optical that an ad¬ 
ditional emitting component is required. They probed different 
spatial scales above the photosphere and found a size variation 
with respect to wavelength in the titanium oxide band ( 0.68 - 
0.82 pm) as well as an asymmetric brightness profile, which 
they attributed to scattering by dust. The indications of asym¬ 
metry in the circumstellar environm ent of L 2 Pup, suggested by 
polarimetry ( iMagalhaes et al.l fl 9861) . as well as the indi catio ns 
of dust scattering by interferometric observations (llreland et alJ 
(2004bh : see also Section ED, might indicate a disk or a com¬ 
panion. 

Very recently, in an effort parallel to our own, iKervella et al.l 
(120141) have presented infrared interferometric data and decon¬ 
volved adaptive-optics images, which show a disk-like structure, 
approximately 700 mas across, transecting the star, opaque at J, 
translucent at H and K, and a thermal dust emitter at L. A faint 
loop extends to the northeast. Kervella et al. attributed the dim¬ 
ming to this disk. L 2 Pup is therefore one of the few AGB stars 
with evidence for asphericities in their circumstellar material, 
and as the closest mass-losing AGB star, is ideally placed for 
resolving the details of the envelope structure. We note that the 
technique of aperture masking used here is different and comple- 
menta ry to the technique of direct adaptive optics imaging used 
bv IKervella et al. ( 2014) . as discussed in detail in Sect. 14.51 

The circumstellar environment of L 2 Pup was observed in 
the near-infrared, using aperture masking techniques on an 8.2 
m telescope (NACO/VLT), to detect any asymmetries within a 
0'.'5 radius. The observations are presented in the next section, 
while results are discussed in Sects. [2 and |4] Conclusions are 
discussed in Sects. 0and[6] 


2. Observations and data reduction 


L 2 Pup was observed on 1 December 2012 (Prog.ID: 090.D- 
0541, PI: D. Klotz) at NACO/VLT dLenzen et a 1 . 12001 
Roussetetak 20031) in sparse apertu re masking mode (SAM; 


Tuthill et al.1120lOt iLacour et al.ll20lH) using the 1 8-hok[] mask 


at 1.64, 2.30 and 3.74 pm. 

Additional observations were obtained on 8 February 2013 
(Prog. ID: 090.D-0677, PI: A. A. Zijlstra) at 2.30 pm using the 
SAM polarimetric mode (SAMPol). In this mode, a Wollaston 
prism positioned beyond the aperture mask splits the incom¬ 
in g beam in to two polarizations (for a detailed description, see 
iNorris et~akll2012l) . With the use of a half-wave plate, the plane 
of polarization of the beam is rotated at 0°, 22.5°, 45°, and 67.5° 
to obtain the polarized intensities at 0°, 45°, 90°, and 135°, that 
is /(>, I 45 , /90 and / 135 . The Stokes vector component, I, of the 
polarized beam can be r ecovered from these int ensities through 
the following relation of lMurakawa et al] ( 2005 1): 


/ - — (/0 + hs + ho + T 135 ), 


( 1 ) 


where T is the transmission of the Wollaston prism (~93% at 
2.30 pm). 

Table Q] contains details of each observing run, such as stel¬ 
lar phases, integration times, number of acquisition frames, and 
seeing conditions. The projected 136 baseline lengths were in 
the range of 0.5 and 7.8 meters, which within only 10 minutes of 
observing time provides an excellent nv-plane coverage (Fig.[D 

2 Although the mask has 18 holes, only 17 holes are used effectively, 
because one is located at the edge of the pupil and is thus omitted in the 
data reduction. 


SAM non-redundant masks allow recovering diffraction-limited 
images with a relative resolution of 20, 30, and 50 mas at 1.64, 
2.30, and 3.74 pm, respectively, with minimal disturbance by at¬ 
mospheric seeing. 

Pupil-tracking mode was enabled to prevent the apertures 
from being covered by the telescope’s spiders. This leads to field 
rotation on the detector, indicated by the rotation of the Mr-plane 
coverage as seen in Fig. [Q All speckle frames were recorded in 
cube mode. To decrease readout time and the volume of recorded 
data, two subframes were chosen on the detector: 1024 x 1028 
pixels for the S13 camera and 512 x 514 pixels for the L27 cam¬ 
era. Different dither positions in all four quadrants of those sub- 
frames were used in a clockwise pattern to minimize readout 
time and detector noise (due to the variation of photon and ther¬ 
mal noise in each quadrant during data acquisition), as well as 
simultaneously recording sky background at an opposing dither 
position. To allow for determining the transfer function and cal¬ 
ibration of the observables, n Pup (K3Ib) was observed immedi¬ 
ately after the science target. 

The data sets were p rocessed with a pipeline developed at th e 
University of Sydney dTuthill et alJ l200ft ILacour et alJl201 fl) . 
The interferograms in each data cube were flat-fielded, dark, 
bias, and sky subtracted. Each speckle frame was Fourier trans¬ 
formed, and the extracted Fourier components (power spectra, 
bispectra) were calibrated with respect to the unresolved target. 
Because of calibration problems, the results from the 1.64 pm 
data are inconclusive. 


3. Fitting the visibilities: indications of an extended 
component 

3.1. Past interferometric observations 

lOchsenbein & Halbwachsl (119821) predicted a diameter for L 2 
Pup in the visual of 13.18 milliarcseconds (mas) based on cal¬ 
culations from the intrinsic color and brightness, llreland et al.l 
(l2004bl) previously studied L 2 Pup with optical interferomet¬ 
ric techniques using aperture masking, although their projected 
baselines were oriented at 180°, that is north-south. They mea¬ 
sured an extended component whose size varies with wavelength 
(Gaussian full-width-half-maximum76±6 mas), and the diame¬ 
ter of a compact source (uniform disk, UD) at 25 + 2 mas (~50% 
of the total flux), which is almost twice the size of the theoretical 
prediction of lOchsenbein & Halbwachsl ( 1982 ). 

VINCI/VLTI data dMeisnerl 120081) show a similar size for 
the K-band emission (2.2 + 0.2 pm\ Pud = 28/)j mas). This 
VINCI diameter varies over less than a month of observations, 
but the visibility is dependent on the different projected base¬ 
lines and on the stellar pulsation phase, that is the duration of the 
VINCI observations spanned over one fifth of the stellar pulsa¬ 
tion period, as shown in Fig. 12 IKervella et af1(l2014l) also found 
evidence for an unresolved component with an upper limit of 
Pud = 17.5 + 1.6 mas, which they interpreted as the central star. 
The resolved components (76 and 28 mas; MAPPIT) would be 
affected by circumstellar emission. 

3.2. This work 

In the SAM power spectra (Fig. |3}, a dependence of the visi¬ 
bility on the baseline position angle (PA.) is shown at similar 
baseline lengths for both wavelengths. The distinction is clearer 
at 3.74 pm, where the difference for visibilities at P.A.< 90°and 
P.A.> 160°is V 2 ~ 0.2 . 
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Fig. 1. SAM/NACO wv-plane coverage with the 18-hole mask for 1 December 2012 at 2.30 /rm (left) and 3.74 /rm (right). Each uv point is 
color-coded to the respective squared visibility value descending gradually from unity (red) to zero (dark blue). 

Table 1. Observation log. 


Mode 

Date 

Phase 

Filter 

d c (A A) [gm] 

# frames 

Camera 0 

tint [sec] 

Seeing* 

Airmass 

SAM 

01 /12/2012 

0.358 

NB_1.64 

1.64 (0.04) 

1200 

S13 

0.5 

1.00 

1.10 

SAM 

01 /12/2012 

0.358 

IB_2.30 

2.30 (0.06) 

600 

S13 

1.0 

0.72 

1.09 

SAM 

01 /12/2012 

0.358 

NB 3.74 

3.74 (0.02) 

800 

L27 

0.5 

1.18 

1.08 

SAMPol 

08/02/2013 

0.855 

IB 2.30 

2.30 (0.06) 

2000 

S13 

0.5 

1.84 

1.07 


Notes. Camera pixel scales are 13.22 and 27.19 mas/pixel for the S13 and L27 cameras, respectively. In units of arcseconds. 


For an initial determination of the overall extent of the de¬ 
tected structures, a simple source morphology was assumed 
(centro-symmetric) and geometric models were directly fitted to 
the Fourier data. The initial sizes were estimated by applying 
a circular Gaussian distribution fit to the azimuthally averaged 
visibilities, and the errors were estimated from the scatter of the 
azimuthally averaged measurements. None of these fits were sat¬ 
isfactory since the dependence of the visibility on the baseline 
P.A. was not taken into account. A further examination of the 
interferometric data was applied within a 20° cone, centered on 
the north-south and east-west directions. Tests were performed 
for both non-averaged and azimuthally averaged squared visi¬ 
bilities. When we applied new Gaussian distribution fits in all 
cases, the radius of the object was larger along the east-west ori¬ 
entation. A similar result is seen in the SAMPol data at 2.30 fim. 
This preferential extension can also be noted in the uv-planes of 
the original data (Fig. [TJ- 

Using the geometric model simulator GEM-FIND (Klotz et al.l 
1201 2 bl) . we found that the shape of the object can be described 
by an ellipsoid, composed of a circular uniform disk represent¬ 
ing the star and an elliptical Gaussian distribution representing 
the circumstellar environment. The best fits to the Fourier data 
are shown in Fig. [4] Simpler centrosymmetric functions are also 
shown for comparison (a uniform disk and a Gaussian distribu¬ 
tion). 

The two-dimensional fitted structures are also overplotted for 
comparison to the NACO images in Fig. [6] The ellipsoid extends 
at approximately 140 mas along the major axis and at a similar 
position angle (~86°) in both wavelengths, while the size of the 


central component (uniform disk) varies from 42 + 0.3 mas at 
2.30 yum to 25 + 0.5 mas at 3.74 gm. 


4. SAM/NACO imaging 

Diffraction-limited maps were reconstructed from the Fourier 
observables using different methods to independently validate 
the results. These included the BiSpectrum Maximum Entropy 
Method (BSMEM: Buscherl 19941) . the MArkov Chain IM- 
ager (MACIM; lire land et al.l 20061) . and the Multi -apertu re Im¬ 
age Reconstruction Algorithm (MiRA; lThiebautll2008l) . Each 
of these methods uses a different algorithm and regularization 
techniques to produce the best-fit map from the interferometric 
observables (i.e., visibilities and closure phases). 


4.1. Image reconstruction with BSMEM 

The reconstruction algorithm BS MEM is based on the maximum 
entropy method (MEM; ISivial 1 19871: ISkilling & Brvanl 1 19841) . 
Based on the geometric model-fitting results (Section [3]), im¬ 
ages were reconstructed from the power spectra and the closure 
phases of the calibrated interferograms using a Gaussian distri¬ 
bution as an initial model. Although the algorithm converged at 
X 2 ~ 1, a small misfit was found at low spatial frequencies (that 
is, baselines shorter than one meter) between the final modeled 
power spectrum and the data at 2.30 yum, suggesting either a mis- 
calibration of the power spectrum at these spatial frequencies or 
that a more complex model was required. 
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Fig. 2. VINCI/VLTI observations of L 2 Pup from 12 October 2001 to 4 
November 2001 with respect to the stellar pulsation phase. The variabil¬ 
ity of the visibility is shown in the top panel, color-coded to indicate its 
dependency on the baseline position angle (25°-70°), while the respec¬ 
tive uniform disk diameters are shown in the bottom panel. A folded 
visual light curve (open circles) is also shown for comparison in the top 
panel. The projected baseline length varies also by less than 1.5 meters 
(14.66-15.98 m) with respect to the physical baseline length (16 m). 


Nevertheless, the image reconstruction produced similar 
structures, both in the SAMPol and SAM data in that wave¬ 
length (Figs.0and[6]). In the SAMPol Stokes I map of February 
2013, two bright clumps appear extended to ~65 mas in the east- 
west direction at levels above 2% of the peak intensity. In the 
SAM maps in particular, a bridge-type structure is more promi¬ 
nent in the west, while a hook-like pattern emerges in the east. 
It is unclear whether the circumstellar environment of L 2 Pup is 
polarized at 2.30 pm. The Stokes I map shown here was pro¬ 
duced by conventional calibration, and this method could not 
provide more information on the source’s polarization. However, 
when usin g a differential calibration method for the SAMPol 
data (INorris et al.!l2012 i). the method is sensitive to subtle po¬ 
larization features. A more detailed analysis of the polarimetric 
data will be presented in a future paper. 

A similar approach was followed for the 3.74 pm data. Here, 
no other misfits to the power spectra were found within the er¬ 
rors. In the maps presented in Fig.[6j we detect an elliptical struc¬ 
ture with sizes 136 mas and 70 mas along the major and minor 
axes. This structure is also aligned in the east-west direction. At 
levels below 5% of the peak intensity, two hook-like patterns are 
seen in the southwest and northeast directions. 



Baseline [m] 




0 18 36 54 72 90 108 126 144 162 180 


Position angle 

Fig. 3. Calibrated V 2 of L 2 Pup at 2.30 pm (top) and 3.74 pm (bottom). 
Each visibility point is color-coded according to the respective baseline 
position angle. 


Extended emission might be resolved out by the instru¬ 
ment, and/or the stellar photosphere is unresolved at a given 
wavelength. We emphasize that the MEM process selects the 
smoothest map that fits the original d ata, based on po sitivity 
and finite-extension of the map dNaravan & Nitvanandalfl98a : 
ISkilling & Brvanll 19841) . and it is thus possible that another so¬ 
lution may have fitted the data better. MEM is also known to 
smooth out fin e stru ctures, for instance embedded point sources 
(iTuthill et al.1120021) . For these reasons, imaging results were re¬ 
produced with other different algorithms, as shown in the fol¬ 
lowing sections. 

4.2. Image reconstruction with MACIM 

To verify the con siste n cy o f the previous reconstructions, we 
used MACIM (llreland et al.l 12006 ). This algorithm uses a reg¬ 
ularization method based on simulated annealing (for a detailed 
description, see llreland et al.ll2006l) . The initial model consisted 
of a uniform disk of 25 mas in diametei[J at 10% of the total flux. 
The results are strikingly similar to those of BSMEM (> 1 % of 
the peak intensity), as shown in Fig. [6] 

4.3. Image reconstruction with MiRA 

MiRA is an algorithm that minimizes a penalty function com¬ 
posed of two terms following a gradient descent method 
(lThiebaut|[2008l) . The two terms are the distance to the data (x 2 

3 Based on the MAPPIT, SAM/NACO and VINCI diameters. 
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Fig. 4. Best geometric model fit (uniform disk + elliptical Gaussian; 
purple) to the SAM/NACO data (gray; as in Fig. [3]i. The fitting results 
are given in the legend. The spread of the values at similar baseline 
lengths is due to the different baseline position angles. 2.30 pm and 
3.74 pm data in the upper and lower panel. Overplotted for comparison 
are the best fits of simpler, centrosymmetric functions, i.e., a uniform 
disk (green; lower line) and a circular Gaussian distribution (blue; mid¬ 
dle dashed-line). The diameters (#ud) and FWFIM of these fits are upper 
panel 223 and 102 mas, and tower panel 180 and 99 mas. 



-150 -100 -50 0 50 100 150 
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Fig. 5. SAMPol Stokes I map at 2.30 pm (February 2013) reconstructed 
with the BSMEM algorithm. Two clumpy structures emerge at ~65 mas 
in the east-west direction at levels above 2% of the peak intensity with 
a possible extension to the northeast. The central star (Lo Pup) is unre¬ 
solved at this wavelength (6 < 30 mas, > 50% peak intensity), but the 
surrounding envelope (6 ~ 50 mas) is non-spherical. Contour levels are 
at 1, 2, 5, 10, 20, 50, and 80% of the peak intensity. North is up and east 
is on the right. 

4.4. Comparing image reconstruction and GEM-FIND fitting 
results 

In the right column of Fig. [6] the image reconstruction results are 
compared to the best-fit models of GEM-FIND. While improve¬ 
ments can be implemented, the models provide reasonably fits 
to the reconstructed images (reduced^ 2 is 17.48 and 6.34 for 
2.30 pm and 3.74 pm, respectively; Fig. [4}: 

Circular uniform disk: the flux ratio of the central star over the 
circumstellar envelope is found to be 0.441 and 0.704 at 2.30 and 
3.74 pm. The respective modeled stellar diameters coincide with 
contour levels at 50 and 80% of the peak intensity in the maps. 

Elliptical Gaussian: the major and minor axes, including the 
position angle of the ellipse, fit the structures in the maps to 
within 5% of the peak intensity (cf. contours in Fig. [6J for both 
wavelengths. 


4.5. Comparison to Kervella etal. 2014 

The observations presented here are taken with the same instru- 
statistics) and the distance to some priors (such as ‘regulariza- ment as was used in lKervella etakl d20H . but with a very dif- 

tion’, which is used to fill the Fourier space between the points of f ere nt technique. Kervella et al. used direct imaging with adap- 

the wv-plane). MiRA implements many different regularizations; tive optics as well as lucky imaging with very short integrations, 

here we used the total variation with edge-preserving smooth- selecting the best images. Deconvolution was used to further im- 

ness with an adopted reg ularizatio n weight p = 10 3 (for a de- prove the resolution. This gives superb resolution but depends on 

tailed description, see lThiebaui2008t). This regularization mini- the accuracy of the point-spread-function (PSF) calibration and 

mizes the total gradients in the image and subsequently smooths 0 n the Strehl ratio. The current data were taken with an aper- 

the flux levels in the image. The BSMEM and MACIM images ture mask. This creates aFizeau inteferometer, giving an exactly 

are similar to MiRA, which succeeds in reproducing the over- _ 

all shape of the structures, including the clumpy and hook-like 4 Although the reduced x 2 statistic of each model fit exceeds unity, 
patterns seen in all reconstructions (Fig. [6]). these were the best fits found. 
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Fig. 6. Comparison of image reconstruction for SAM maps at 2.30 pm (top) and 3.74 pm (bottom) from December 2012. Each column, from 
left to right, represents image reconstructions with BSMEM, MACIM, and MiRA. Contour levels and orientation are as in Fig. [5] The resolution 
elements of SAM for the respective wavelengths are shown in the lower left corner of the second column. GEM-FIND model fits are overplotted in 
the third column for comparison: elliptical Gaussian (FWHM; orange) and circular uniform disk (white). The extent and size of the fitted models is 
similar to those of the mapped structures down to 5% of the peak intensity, while >70% of the emission originates from the central star. At 2.30 pm 
two clumpy patterns emerge in the east-west direction with a possible hook-like pattern rising in the north-east. The former are unresolved at 
3.74 pm, and the shape now resembles a bar with two hook-like patterns emerging in the north-east and south-west directions. 


known PSF that is independent of the seeing (the seeing does 
affect the signal-to-noise), but with reduced throughput because 
of the mask geometry. 


Comparing the two techniques, aperture masking has higher 
angular resolution: its PSF is somewhat smaller than the diffrac¬ 
tion limit of a single mirror and deconvolves well because it is 
mathematically known. Both techniques require accurate cali¬ 
bration (PSF for direct imaging and visibility amplitude for aper¬ 
ture masking). Aperture masking filters out extended emission, 
while adaptive optics leaves extended pedestals in the PSF. The 
techniques are complementary in view of their very different 
strengths. 


Comparing our Fig. [5] with Fig. 4 of iKervella et al.l (l2014h 
illustrates this. The current data have a slightly better resolution 
and show smaller details, but the extended disk is better seen in 
Kervella et al. Allowing for this, the agreement is very good. The 
aperture masking K -band image shows two clumps either side of 
the star, indicating a better resolution than Kervella et al., but the 
extended disk emission is suppressed. The L-band image, with 
its shorter uv-spacings, shows the disk better. The overall struc¬ 
tures agree very well, validating the images of iKervella et al.l 
( 2014 ) and supporting the reliability of the image reconstruction. 


5. Nature of the detected structure 

The NACO observations are focused at 2.30 pm, where the de¬ 
tection of carbon monoxide emission^ and scattered dust emis¬ 
sion from the circumstellar envelope is expected, and at 3.74 pm, 
where continuum emission from the dusty environment probably 
is more prominent. The detected ellipsoidal structure is com¬ 
mon in both wavelengths. This was confirmed by convolving 
the 2.30 pm maps with the resolution beam of SAM/NACO at 
3.74 pm. The overall structure may be interpreted either as an 
elongated outflow or a toroidal structure, depending on the in¬ 
clination of the object with respect to the line of sight and the 
stellar rotation axis. The interpretation of the hook-like features 
depends on the inclination of the object as well. A few possible 
scenarios on the origin of the structure detected around L 2 Pup 
by SAM/NACO are presented below. 

5.1. Circumstellar shell hypothesis 

The structure detected at 2.30 pm is aspherical and extends on 
both sides to 65 mas on average (above 3% of the peak flux east- 
west), which in physical units translates in 4.2 au. The structure 
becomes an ellipsoid at the pseudo-continuum (3.74 pm), with 


5 To our knowledge, there are currently no spectroscopic observations 
of the CO first-overtone band available in the literature. 
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approximate sizes 140 mas wide in the east-west direction and 
70 mas w ide in the north-south direction (Fig. [6]). 

IBedding et al.l (12002}) compared the IRTS mid-infrared spec¬ 
trum (recorded in 1995, pre-dimming) to a model of a silicate- 
rich, detached shclQ The inner and outer radii of that dust shell 
were 7.2 x 10 14 and 1.3 x 10 15 cm (or 48.1 and 86.9 au). At 
a distance of 64 pc, these dimensions translate into 0775 and 
1735, and are placed within the extended envelope detected at 
11.7 yum (178x178: lJura etari|2002 ). With_a known expansion 
veloc iuQ of approximately 2.5 km s 1 dKerschbaum & Olofssonl 
ll999l) . the outer edge of the shell would have reached that dis¬ 
tance within 165 years. Nevertheless, such a shell is much larger 
than the structures detected by S AM and also lies outside the 
region imaged bv lKervella et al.l ( 20141 . 


There have been previous detection s of molecular shells in 


the extended atmospheres of AGB stars ( 

Le Bouquin et al. 

2009; 

Perrin et al. 2004 Eisner et all 20071; 

Thompson 

et al. 

2002; 

Ohnaka et al. 20ll IWittkowski et al. 2007, 201 It 

Ireland et al. 

2004a). The sizes of these shells have been estimated through 


star+molecular shell models (MOLspheres) or dynamic atmo¬ 
sphere models and are spherically symmetric. The observations 
presented here showed that the case of L 2 Pup is different be¬ 
cause no spherically-symmetric shell was detected. 


5.2. Toroidal structure hypothesis 


The molecular emission of L 2 Pup has been exten 
sively monitored in the past fifty years; it showed the 
variab ility of the outflow velocity 

IKerschbau m & Olo fs sonl 11999} 

2002} [Lenineet al.l 1976} Winters et al 


(IKerschbaum et af 


1990} lHaikala et al. 


Winters et al. 

20021 


2003 


iMenten & Melnic' 


Haikala 


1991; 


Gomez Balboa & Lepine 


11994 

S- 

atchelor 


Lepine et al. 1978j 


ck 

I 


etak[ 197l[| Knowles & Batchgloi[ 1978^ Bahsteretall 


Caswel 


.1977; Schoieretaf 12004 iGonzalez Delgado et al.l 20031) . 
Gonzalez Delgado et al.l (120031) reported that the interior of the 


circumstellar envelope of L 2 Pup is probably complex, based 
on their SiO obser vatio ns. However, the radio interferometric 
measurements by ISchoier et alJ (120041) of “thermal SiO” were 
unable to probe the direction of these flows, and a circularly 
symmetric envelope was assum ed. 

Recent observations (iMcIntosh & Indermuehlell2013l) of the 
variation in the positioning of the velocity centroids of SiO 
maser emission over a period of four years, indicate an asym¬ 
metric nature of the gaseous environment of L 2 Pup which 
may be associated with a rotating edge-on disk and bipolar 
ejecta (see Sect. 15.31 ). SiO maser emission probes the mate- 
ri al muc h closer to the star than the “thermal” SiO emission. 
IKerschbaum & Olofssonl ( 1999 ) suggested that the line profile 
seen in the carbon monoxide spectra is caused by a gaseous disk 
around the star, s i milar to the case of SV Psc dKlotz et al.l2012ah . 
iBuiarrabal et al.l (12013l) showed that similar CO line profiles can 
be found in post-AGB stars, and they assigned these to rotation 
in molecular disks. In fact, three of these objects - the Red Rect¬ 
angle, AC Her, and HR 4049 - have slow expansion velocities, 
bipolar ejecta, and circumbinary dusty disks. 

The rotational axis of L 2 Pup is unknown, which allows for 
a range of scenarii on the nature of the detected structures. How- 


6 Nevertheless, they advised to use the best-fitting parameters of their 
model with caution. 

7 This is the outflow velocity of gas, not dust. These are proba¬ 
bly only loosely coupled in a star with such a low mass-loss rate 
(~ 3 x lO^Moyr 1 ; ISchoier et al.ll2004) . 


ever, the images of lKervella et al.l d2014l) are strongly suggestive 
of a rotating east-west disk, which would indicate a north-south 
rotation axis. 

Considering that the bilateral clumps are observed farther 
out than the extended circumstellar emission (^30 mas) in all 
2.30 yum maps, but are absent from the 3.74 pm mapfQ it is 
possible that the edge-on structure detected with SAM/NACO 
is in fact a torus of gas interwoven with dust. The origin of the 
hook-like patterns is quite perplexing. These emerge at levels 
<5% of the peak intensity and extend symmetrically southwest 
and northeast in all the SAM maps. Assuming that the inclina¬ 
tion of L 2 Pup in our line of sight is perpendicular to the north- 
south axis, that is 90°, the possibility of an inclined, edge-on 
spiral or helical st ructure, similar to those ob served in AGB stars 
(e.g., AFGL3068: lMauron & Hug gins||2006|) and planetary neb¬ 
ulae (e.g., PNG356.8+03.3; ISahai et al.ll201 ll) . should not be ex¬ 
cluded. 


5.3. Stellar and torus parameters 

The SiO maser spectra of lMcIntosh & Indermuehlel(l20L3l) show 
an unusual profile, with two peaks separated by 20 km s’ 1 , in ad¬ 
dition to the more usual narrow peak at the systemic velocity. A 
rotating gaseous disk would fit this profile, where the separated 
peaks show the approaching and receding side of the disk. The 
central velocity component may trace the stellar wind might be 
emission from the disk in front of the star, where the maser is 
amplifying the stellar emission. This model gives a rotational ve¬ 
locity of the disk of v sin i = 10 km s’ 1 , where i is the inclination 
with respect to the plane of the sky. 

The two clumps seen on either side of the star are consistent 
with the expected structure seen at the inner radius of the disk. 
This gives an inner radius of 4.2 au (Kervella et al. derived 6 au 
from their models). 

Assuming Keplerian velocities, this gives a mass of the star 
of M sin /’ - 0 . 7 M f . ; . The optical obscuration events, and the im¬ 
ages of lKervella et al.l (120141) . indicate an inclination i ~ 90°, 
although given the size of the star, a significant range of i is 
possible. However, using a maximum range of i = 60 deg and 
r - 6 au still implies a stellar mass M < 1 M e . 

This mass is lower than that derived by Kervella et al. (2 M e ). 
They located the star at the beginning of the AGB (meant is 
the start of the TP-AGB) where the mass is still close to the 
initial mass. However, the short period of the star of 140 days 
together with a location 0.5 mag below the M ira PL relation 
(IBedding & Ziilstral 11998} IBedding et al.l 120021) suggest a low 
initial mass, as found here. 


5.4. Binarity hypothesis 

The possible signature of a binary companion was investigated 
in the SAM/NACO data, a lthough the bin ary na ture of L 2 Pup 
was previously doubted bv ljorisse net al.1 ( l2009t>. Using the ge¬ 
ometric model-fitting tool LITprdl (iTallon-Bosc et alil2008l) . an 
additional, unresolved component, with a flux ratio ranging from 
1/100 to 1/500, was included in the model presented in Sect. 0 
and was positioned randomly within the NACO field of view. All 
attempts failed to converge (x 2 >> 1), suggesting that no signal 
from a binary companion was found in the data. 


8 Due to the lower resolution at the wavelength. 

9 LITpro software available at http://www.jmmc.fr/litpro 
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6. Conclusions 


The SAM and SAMPol diffraction-limited imaging revealed the 
complex circumstellar nature of the envelope of L 2 Pup in both 
gaseous emission and dust continuum. The maps at 2.30 pm 
show an asymmetric structure extending to ~65 mas in the west 
at levels above 5% of the peak intensity and a hook-like exten¬ 
sion (Figs. [5]and[6]). It resembles a linear structure from the star, 
whose extended envelope is resolved at ~30 mas, although the 
stellar photosphere should be smaller than that. A fainter, hook¬ 
like feature emerges in the northeastern direction. The overall 
shape of the structure is similar at 3.74 pm, where it resembles 
an ellipsoid with sizes 140 mas and 70 mas along the major and 
minor axis, respectively. This warm dust component extends far¬ 
ther out into a fainter, helical feature. 

Both structures are oriented in the east-west axis, agreeing 
with iKerve 11a et all (120141), but different from previous aperture- 
masking observations: Ireland et al.l (l2004bl) positioned their 
mask in a north-south orientation (PA. = 180°) and failed to de¬ 
tect the equatorial structure. 

This ellipsoidal structure, found in pseudo-continuum 
(3.74 pm), is the less resolved counterpart of the more detailed, 
bilateral structure (2.30 pm), since the dimensions of both are 
similar. The structure is co ntained within the extended silicate- 
rich dust zone dBedding et al.l 120021) . As shown in the near- 
infrared maps, it appears that the structure detected in L 2 Pup 
is seen nearly edge on. We interpret the structure as a torus with 
inner radius 4.2 au around L 2 Pup. This toru s is probably the 
cause of the long-lived dimming event (iBedding et al.ll2002l) . It 
might be that a more extended dust reservoir is locate d between 
the torus and the extended envelope djura et al.ll2002 ). although 
it is currently outside the field of view of SAM/NACO. 

SiO maser profiles indicate that the torus may be rotating at 
10 km s -1 . Using the inner radius of 4.2 au, we derive a stellar 
mass of 0.7 M e with an upper limit of 1 M Q . This low mass is 
consistent with the short period of L 2 Pup. 
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Appendix A: Visual light curve 

Figure IaTI shows the visual light curve of L 2 Pup from 1985 to 
2013 (data from AAVSO), where the observing periods of MAP- 
PIT, VINCI, and NACO are clearly marked. An average light 
curve is overplotted to guide the eye. 
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Fig. A.l. Visual light curve of L 2 Pup from 3 January 1985 to 6 July 2013 (data from AAVSO). The dimming event started in 1994, and the stellar luminosity in the V band has not increased to 
pre-1994 values ever since. The observing periods for MAPPIT, VINCI, and NACO are marked in magenta, green, and blue lines, while the selected Julian date for zero pulsation phase, selected 
on JD=2 446 322, is drawn as a red box. An average light curve (gray) is overplotted to guide the eye. The inset shows a magnified version of the fitted light curve over two pulsation cycles, during 
and right after the MAPPIT observations. 
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